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We present experimental results on two-qubit Rydberg blockade quantum gates and entanglement 
in a two-dimensional qubit array. Without post selection against atom loss we achieve a Bell state 
fidelity of 0.73±0.05, the highest value reported to date. The experiments are performed in an array 
of single Cs atom qubits with a site to site spacing of 3.8 ^m. Using the standard protocol for a 
Rydberg blockade Cz gate together with single qubit operations we create Bell states and measure 
their fidelity using parity oscillations. We analyze the role of AC Stark shifts that occur when using 
two-photon Rydberg excitation and show how to tune experimental conditions for optimal gate 
fidelity. 
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I. INTRODUCTION 

Qubits encoded in liyperfine states of neutral atoms 
are a promising approach for scalable implementations 
of quantum information processing 0. We are develop¬ 
ing an atomic qubit array for quantum logic experiments. 
The array consists of qubits encoded in Cs atom hyperfine 
states. Single qubit gate operations are performed using 
either microwave fields for global operations on the array, 
or focused light fields for control of individual qubits [2|. 
Two-qubit entangling gates are based on Rydberg block¬ 
ade interactions [3j. Qubit initialization is performed with 
optical pumping and qubit readout is based on imaging 
of resonance fluorescence Q. 

Provided sufficiently high gate fidelities can be 
achieved the neutral atom approach provides a scalable 
path towards large qubit numbers. The qubit density in 
our recent 2D implementations 0, 0 is approximately one 
qubit per 14 /im 2 with a loading fraction of 60%. This 
translates into an effective area per qubit of 24 p, m 2 . The 
area needed for a large number of qubits, say 10 6 , would 
be a modest 0.24 cm 2 . Although there are numerous 
engineering challenges associated with scaling to such a 
large number of qubits there is no fundamental reason 
why this could not be achieved. 

At the present time the largest impediment to scaling 
is that the demonstrated gate fidelities have not reached 
the level where fault-tolerant coding and error correc¬ 
tion architectures are viable [gj]. Single qubit gate oper¬ 
ations have reached better than 0.99 fidelity with single 
site control 0 and it is reasonable to anticipate further 
improvement using composite pulse sequences (tJ. The 
fidelity achieved to date for entangling gates is less sat¬ 
isfactory. Three research groups have demonstrated en- 
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tanglement of neutral atom qubits using Rydberg inter¬ 
actions. The results have been characterized in terms of 
Bell state or entanglement fidelities with reported val¬ 
ues of 0. 580, 0.710, 0.75 pH, and 0.81 (Til allowing for 
post selection to correct for atom loss during the gate se¬ 
quence. Reported fidelity results without post selection, 
which is preferable for quantum computing applications, 
are 0.580 and 0.60fllj. 

The experimental entanglement fidelity results re¬ 
ported to date lag far behind theoretical analyses which 
predict gate fidelities > 0.99 [HJ in a room tempera¬ 
ture apparatus and > 0.9999 for circular Rydberg states 
at cryogenic temperatures [l3j. It has therefore been 
an open question as to whether or not the separation 
between experimental and theoretical results is due to 
purely technical errors, or derives from some unaccounted 
for aspect of the atomic physics. We demonstrate here 
that previous analysis has not fully accounted for AC 
Stark shifts that occur in two-photon excitation of Ry¬ 
dberg states. We clarify the impact of the Stark shifts 
on the effective gate matrix, and show how to minimize 
sensitivity to imperfectly controlled experimental param¬ 
eters. We then demonstrate improved two-qubit entan¬ 
glement with fidelity of 0.79 ± 0.05 allowing for post se¬ 
lection and 0.73 ± 0.05 without post selection. Although 
still below what is needed for scalability we anticipate 
that further improvement will be possible in the future. 

The rest of the paper is organized as follows. In Sec. 
El we recall the Rydberg blockade Cz protocol and an¬ 
alyze the impact of AC Stark shifts on the gate. We 
then proceed to show how to compensate for the Stark 
shifts to obtain an ideal gate matrix. In Sections IIIIA . 
B we describe the experimental setup and how to mea¬ 
sure relevant parameters using Ramsey interference and 
Rabi oscillation experiments. In Sections IIIIC .D we de¬ 
scribe CNOT and Bell state experiments, followed by a 
concluding section IIVI with an outlook on future devel¬ 
opments. 
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FIG. 1. (color online) Level structure (left) and pulse se¬ 
quence (right) for two-qubit Rydberg blockade Cz gate. The 
qubits are encoded in ground hyperfine states |0), |1) while |r) 
is a high lying Rydberg state. 
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FIG. 2. (color online) Level diagram (a) and ac Stark shifts 
(b) for two-photon Rydberg excitation. 


The standard protocol for creating a controlled phase 
gate via Rydberg blockade uses a three pulse sequence to 
implement 


/I 0 0 0 \ 

0-10 01 
0 0-10 
\0 0 0 -1 / 


(1) 


The sequence shown in Fig. |T] uses a 7r pulse on the 
control qubit, 27 t on the target, and n on the controlQ. 
In the ideal situation of perfect blockade and negligible 
ratio of excitation Rabi frequency Q to qubit frequency 
splitting w q we obtain the gate matrix ©■ While sev¬ 
eral experiments have used single photon excitation of al¬ 
kali atom Rydberg states [3|,[HI, all the Rydberg based 
quantum gate experiments except [ll| have used a more 
complicated two-photon excitation method. The primary 
reason for doing so has been to avoid the need for high 
power at the short wavelengths of one-photon excitation 
(297 nm in Rb and 321 nrn in Cs). As we proceed to show, 
the use of a two-photon drive changes the gate matrix so 
that even under ideal conditions we do not obtain Eq. 
©• 
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with cj> shorthand for the phases {<(>oo> (f>oi,<fiio, <j> n}. The 
C z ^ operator can only create entanglement when 0oo — 
<A)i ~~ 4>w + <^n 7 ^ 27m with n integer. It is therefore 
essential to correctly control the gate phases. 


A. AC Stark shifts 

The first contributions to the gate phases come from 
ac Stark shifts that arise due to the use of two-photon 


excitation as shown in Fig. [2] We divide the Stark shifts 
into resonant, and non-resonant contributions. Each of 
the qubit levels |0), |1) acquire resonant and non-resonant 
Stark shifts as does the Rydberg level. All shifts are listed 
in Table H] 

There are several things to note about the expressions 
given in Table Q] The resonant shifts are the standard 
expressions valid for the situation where |Ai| y e with 
7 e the radiative linewidth of the intermediate level. We 
assume two-photon resonance between |1) and |r) so that 
A 2 = — Ai and A = Ai-FA 2 =0. The field amplitudes 
and Rabi frequencies are related by flj = dj£j/h with dj 
the relevant transition dipole matrix element, Ej the elec¬ 
tric field amplitude, and the two-photon Rabi frequency 
LIr = ! 2 ifl 2 /( 2 Ai). In the approximation that the hy¬ 
perfine splitting of |e) is small compared to the detuning 
Ai the expressions given are valid. When this is not the 
case the expressions for the Stark shifts as well as the 
relation between Q and the one-photon Rabi frequencies 
have to be modified. The full expressions for the specific 
case of Cs atoms excited via 6 si/ 2 —> 7pi/2 -A ns i/ 2 , as 
in our experiments, are given in Appendix lAl 

The non-resonant polarizabilities a g j can be calculated 
using a sum over states approach. Since we explicitly 
account for the resonant contributions to the ac Stark 
shift, the polarizabilities in this paper are defined with 
the resonant transitions excluded from the sum. The 
non-resonant Rydberg polarizability is given by the ex¬ 
pression 
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where e and m e are the electron charge and mass, re¬ 
spectively, and ujj is the frequency of field j. We will 
assume that the excitation beams are large compared to 
the size of the Rydberg wavefunction and ignore correc¬ 
tions to the Rydberg shift arising from finite beam size 
effects [3]. 
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TABLE I. Stark shifts contributing to the gate phases. Here £j is the field amplitude of field j, f lj is the Rabi frequency, Ai 
is the detuning from the intermediate level, and a g (n)j is the non-resonant polarizability of the ground(Rydberg) levels at the 
frequency u)j of field j. Superscripts r and nr label resonant and nonresonant Stark shifts respectively. Full expressions that 
account for the hyperfine structure of the intermediate level are given in Appendix lAl 


description 

shift on |0) 

shift on |1) 

shift on |r) 

resonant shift from £\ 

non-resonant shift from £\ 

resonant shift from £2 

non-resonant shift from £2 

A r — 1 l^ll 2 

A 01 - 4 Al _„ q 

Ag = -&a gl \£ 1 \ a 

Ag2 = -^a g 2|£ 2 | 2 

A r _ 1 1^1 I 2 

A H - 4 ~KT 

A nr 

A gl 

A nr 

A g2 

Ari = —3^Q!Rl|fl| 2 

Ar _ 1 | O 2 | 2 

^R2 — J—aT 
Ar 2 — —^Q'R2|f’ 2 | 2 


B. Phase shift from a 2n rotation 


In the case of a two level system driven by a single field 
the ground state accumulates a shift of e I7r = — 1 during a 
resonant 2tt pulse. The situation is more complicated for 
the three level system driven by two fields. The laser fre¬ 
quencies are tuned to give full population transfer from 
the ground state to the excited state. This implies that 
the detuning compensates the Stark shifts from each of 
the excitation beams. In the case where there is two- 
photon resonance, or near resonance, and Ai/ 7 e is large 
to minimize spontaneous emission the resonance condi¬ 
tion is 


A * (A R2 + A& + AS) - (A r n + A£ + A”) = 0. 

(4) 


This includes the resonant Stark shift, which can be can¬ 
celled by setting |fli| = |1 (when the intermediate level 
hyperfine structure is negligible), and the non-resonant 
shifts on the ground (Rydberg) states, A^ R ^ 2 ,. 

In the far detuned limit when the hyperfine structure 
of the intermediate state can be neglected the phase accu¬ 
mulated by the ground state after a resonant 2n Rydberg 
pulse is 




H-i 


4>R = 7T 

1 - 

£2 2 

sign(Ai) 


(A" 1+ A^)2t w , (5) 


with t^ir = 2tn = 27r/1 S~2 r |. Equation J5]) is readily de¬ 
rived using the Schrodinger equation for a three-level 
ladder configuration, and adiabatically eliminating the 
intermediate state in the limit of large detuning. 

We can express the Rabi frequency ratio as 
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The 27 t pulse phase can therefore be written as 


0 R = 7r-(A r 11 + A^+Ag)2t 7r . (6) 

This way of writing the phase has a clear physical in¬ 
terpretation. The factor of 7r is the quantum phase ac¬ 
cumulation from rotating the effective two-level system, 


which is analogous to a spin 1/2, through 2w. The second 
term is the Stark phase accumulated by the ground state 
over a time of 2t^. 

It may be surprising that the usual picture of a 2-7 t Rabi 
pulse imparting a 7r phase shift is only valid for a single 
photon transition. For a two-photon drive the phase shift 
can take on any possible value although values of the Rabi 
frequencies for which </ R = ir can always be found. In 
particular when the ac Stark shifts on the ground state 
are fully compensated, i.e. A| x + A/) + A“£ = 0, then 
</>r = 7r, and we recover the one-photon transition result. 


C. Gate phases 

The C z z operator is constructed by considering how 
the shifts described in the previous sections affect the 
computational basis states, | ct) = {|00), |01), 110), |11)}. 

For the |00) state the excitation beams are detuned by 
w q for both qubits so both remain in the 100) state. From 
Table |T1 the shifts on |0) include the resonant £\ shift and 
the non-resonant £\ and £2 shifts. The shift on |0) for 
the control(target) qubits is 

</hf,c(t) = ~ (^01,c(t) + ^gl,c(t) + ^g2,c(t)) 2f7r, c (t)j (7) 
so 

<A)0 = </> hf,c + <t> hf,t- (8) 

Here we have allowed for a possible variation in param¬ 
eters at the control and target qubit sites so that 
</>hf ,t need not be equal . 

For the |01) state, the control experiences the off- 
resonant Stark shift, </>hf,o while the target picks up a 
phase shift from the resonant 27 t Rydberg pulse, </> R from 
Eq. (j5j), 

4>01 = 0hf,c + (9) 

The 110) state is different than the |01) state because 
the control atom is held in the Rydberg state for a time 
2t g ap + 2t 7T x. The phase accumulated during this time is 
due to the ground-Rydberg differential Stark shift and is 
given by 

0 gap = — [(a R2)C + a r1)C + a R 2 iC ) 

- (A r lljC + A” c + A” c )] 2 (t gap + 1 ^). (10) 
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Here f gap is the extra time in between pulses which, ex¬ 
perimentally, is the minimum time it takes to switch the 
laser beams between control and target sites. In total the 
shift on the 110) state is 

010 = 0R,c “b 0gap 0hf,t- (H) 

Finally the 111) state experiences an additional shift 
due to the blockade, 0 b, for a time 2t v t which includes 
the resonant Stark shift from E\ and the non-resonant 
shifts on the |1) state of the target atom 


To proceed we note that we can always ensure 0r = mr 
by choosing the correct value for the Rabi frequency ra- 

tio q = ^ . To see this let = a|fii| 2 , A™[ = £>|T2!| 2 , 

A ^2 = c|fi 2 | 2 , r = 7rd/|H 1 H 2 | with a,b,c,d real con¬ 
stants that depend on the detuning from the intermedi¬ 
ate level and atomic structure parameters. If we then 
choose q such that 


1 — n = 


(a + b)d 


( cd)q 


0b — (A^i t + AgJ t + Ag^t) 2t Wi t + 0 bl> 

= —7r + 0 Rit + 0 bl- 

The last term is a small blockade leakage phase)!} 0bl = 
7tHr/(2B). The total phase accumulation on the 111) 
state during the gate sequence is thus 

011 = — 7F + 0 R,c + 0R,t + '/’gap + 0BL- (12) 

Equations (l8lfl2l) fully determine the phases of the C z ^ 
operator which we summarize here for convenience 


we obtain 0r = nir. We then set f gap such that 0 gap = 
2n'n. Solutions occur at 

^■gap T ^7r 

ti'tt 

+ A ri + A R 2 ) - (A r n + + Ag) 

Setting 0 r, 0 gap to be multiples of 7r as described 
above, the gate phases modulo 27r, for n odd and n' even 
are 


0 00 = 0 hf,c + 0hf,t, 

001 = 0 hf,c + 0R,t, 

010 = 0R,c T 0gap “b 0hf,t> 

011 = “71 + 0R,c + 0R,t + 0gap + 0BL ■ 

The gate phases are not completely independent since 

001 + 010 — 000 — 011 = 7T — 0BL ■ 

In the limit where the blockade leakage phase 0bl is 
small, which will be the case for parameters which yield 
high fidelity entanglement, there is a fixed constraint be¬ 
tween the phases. A global multiplicative phase factor 
is irrelevant, leaving two free phases. As we discuss in 
the following section a correct choice of two parameters 
is sufficient to fix the two phases and obtain an ideal gate 
operation. 


D. Setting parameters to recover an ideal Cz gate 

In general the C z ^ operator does not necessarily cre¬ 
ate entanglement, and does not directly create Bell states 
with standard phases. In this section we show that it is 
in principle possible to choose parameters such that an 
ideal Cz operator is implemented by the pulse sequence 
of Fig. [I] For simplicity we assume that Hi, H 2 are the 
same for both control and target atoms so that we can 
drop the c, t subscripts on the gate phases. We will also 
neglect 0bl since it is a small error for typical experi¬ 
mental parameters of B > Hr. Alternatively the 0bl 
phase can be cancelled using a slightly modified pulse se¬ 
quence, which does not change the other gate phases, as 
described in Fig. 3 of Ref. 1121 . 


0oo =20hf, 

0Ol=0hf+7T, 

010 =0hf + 7T, 

011 =7T. 

For n even and n' odd we get 

0oo =20hf, 

00! =0hf, (M) 

010 =0hf) 

011 =7T. 

The 0hf phases can be corrected by applying global R z ( 9 ) 
rotations with 9 = — 0hf which recovers the ideal Cz of 
Eq. (QJ. We emphasize that an ideal Cz gate is recov¬ 
ered apart from errors due to spontaneous emission from 
the intermediate and Rydberg levels, finite temperature 
Doppler and position fluctuations, and finite blockade 
strength. Those errors have been quantitatively studied 
in previous work[l2j, but without the constraints implied 
by the parameter choices presented here. We defer a re¬ 
examination of the theoretically achievable gate fidelity 
in a real atom to future work. 


E. Setting parameters to recover a Cx gate 

In the experiments described below in Sec. IHIl we have 
not implemented the parameter settings needed for an 
ideal Cz gate. Nevertheless we can still create a modi¬ 
fied CNOT = Cx gate and create Bell states, albeit not 
with the standard phases. The standard Cx gate in the 
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computational basis {00,01,10,11} is 


Cx 


/l 0 0 o\ 
0 10 0 
0 0 0 1 
\0 0 1 Oy 


An equivalent operator, but with the X operation condi¬ 
tioned on the control qubit being in state |0), is 


C x = (X® I)C x (X ® I) 


/o 1 0 o\ 
10 0 0 
0 0 10 
\0 0 0 ly 


The usual method to transform the Cz gate into a 
Cx gate is to apply Hadamard gates on the target qubit 
before and after the Cz operation. We generalize this 
to 7 t/ 2 rotations, with a relative phase 9 to create the 
operator 


C x ,m = U/2, 9 ) Cz ,} R *{ tt/2, 0) 

where /? c ( t ) (£, 9) is a ground state rotation on the con- 
trol(target) by an angle £, with phase 9. For arbitrary 0 
and 9 we find 


C*,#) 


fa b 0 o\ 

c d 0 0 
0 0 e / 
y0 0 g h) 


with eight nonzero elements a — h. However, for specific 
values of 9 we get a C x ^{9) operator proportional to Cx 
or Cx, with only four nonzero elements of unit modulus, 
but with different phase factors. We define a matrix over¬ 
lap as 0(C X , C X j(0)) = \ E Cx|Cx,^)| 2 . It then 

elements 

follows that 

0(Cx,C x ^(cpu - 0i O )) = 1, 

0{C X , C x ^{<j>oi — 0oo)) = 1. 

Note that the overlap is not a gate fidelity and does not 
contain any phase information. The difference between 
the phases used to prepare the two C x gates is indepen¬ 
dent of parameters as expected, 


(011 — 0io) _ (^01 — 000 ) = — 0R,t + 0B 

= — 7T. 


(16) 


We then use the C x ^(9) gate to create Bell like states, 
but with nonstandard phases. To do so we apply the 
sequence 

U{6,j)=C X j{0)R c W2,0). (17) 


Operating with U on the product state |00) creates a 
maximally entangled state when 9 = 0n—0io- In general 
this will not be one of the Bell states due to the presence 
of different phase factors on the two components of the 
state vector. If desired we can recover a standard Bell 
state with additional one qubit rotations. This is the 
approach we demonstrate in the next section. 

To determine the entanglement fidelity of the state 
created with U we measure the populations and the 
two-qubit coherence terms which can be done by the 
method of parity oscillations [13 • Defining the parity sig¬ 
nal P = Poo + Pu ~ Poi — Pi 0 , where P-ij are the diagonal 
terms of the two qubit density matrix, a 7r/2 rotation at 
an angle 9 on both qubits transforms the parity signal to 

P' = 2Re[C 2 ] -2|Ci|cos(26» + 0). (18) 

Here C\ = \Ci\e 1 ^ is the coherence between states |00) 
and 111) and C 2 is the coherence between states |01) 
and 110). The coherence C 1 can then be extracted 
from the parity oscillation measurements. The entan¬ 
glement fidelity of states close to the Bell state 

is F = Bw+fn _|_ | C ' : |. Values of F > 0.5 are a sufficient 
condition for the presence of entanglement fl8[ . 


III. EXPERIMENT 
A. Set-Up 

Cs atoms are loaded into a 49 site array of blue-detuned 
dipole traps formed by 64 tightly focused, weakly over¬ 
lapping 780 nm beams as described in Refs. The 

7x7 site array has a 3.8 pm site to site spacing. The 
qubits are encoded in the hyperfine clock states with 
|0> = | 6 si/ 2 ,/ = 3,771/ = 0) and |1) = | 6 si/ 2 , / = 
4, rrif = 0). Single qubit rotations are performed glob¬ 
ally using a 9.2 GHz microwave field. Single site rotations 
use a combination of the microwave field and a tightly fo¬ 
cused beam detuned by -14 GHz from the 7pi/z line to 
induce a differential Stark shift equal to ~ 40 kHz on a 
single site. In contrast to the approach of Ref. 0 where 
the microwave frequency was detuned from w q and the 
optical Stark shift provided a local resonance condition, 
here the microwave frequency is set to w q and the Stark 
shift is used to tune the site where no rotation is desired 
out of resonance. 

Rydberg excitations are driven using a two-photon 
transition through an intermediate lpi /2 level. The two 
wavelengths are 459 nm and 1038 nm. The lasers are 
locked to high finesse ULE resonators for long term fre¬ 
quency stability and short term linewidths < 500 Hz on 
few ps timescales. We measure Ai, the detuning from 
the center of mass of 7pi/2, by maximizing the light scat¬ 
tering from the |7pi/2,/ = 4) level, and then detuning 
the light by a known amount. The hyperfine constant of 
7pi /2 is reported in Table Hill The duration of the ex¬ 
citation pulses is controlled by acousto-optic modulators 
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(AOMs). The pulses were of square shape as indicated 
in Fig. [T] with typical rise and fall times of 50 ns. For 
each Rydberg pulse the 1038 nm light was left on for ap¬ 
proximately 50 ns longer than the 459 nm light so that 
the precise value of the pulse duration was controlled by 
the 459 nm light. This leads to some additional ground- 
Rydberg AC Stark shifts from the 1038 nm light that are 
not accounted for in the analysis of Sec. nm 

Each beam is sent through separate fibers to a 2D 
beam scanner which is created using two crossed AOMs 
which allows for 2D positioning across the array 0 - 
The two counter-propagating beams are focused to waists 
(1/e 2 intensity radii) of 3.0 and 3.7 ^m for the 459 and 
1038 nm light, and aligned onto a single site. The site 
to site switching time is ~ 0.5 /zs. The two photons are 
<r+, (7— polarized with respect to the quantization axis 
z which is perpendicular to the plane of the qubit ar¬ 
ray. A 0.15 mT bias magnetic field is applied along z. 
With these polarizations and choice of intermediate level 
we excite a Rydberg \ns\/ 2 , m 7 = —1/2) fine structure 
state. As shown in Appendix [B] there is negligible ex¬ 
citation of the rrij = +1/2 Zeeman state. All data re¬ 
ported here are for the | 82 si/ 2 ,mj = — 1 / 2 ) state with 
Ai = 2tt x 0.83 GHz. The optical trap array is turned off 
for the few /j,s duration of the Rydberg gates. Turning 
the traps back on after Rydberg excitation leads to pho¬ 
toionization or mechanical loss of the atoms before they 
decay to the ground state. In this way trap loss is used 
to measure the Rydberg excitation probability. 

As we discuss in the next section the detuning Ai and 
the beam powers are chosen to minimize the differential 
ac Stark shift between ground and Rydberg states which 
can be expressed as 

A d iff, g R = (A ^ 2 + A& + A&) - (A r n + A“i + Ag) . 

Minimizing this shift is advantageous as it reduces sen¬ 
sitivity to intensity fluctuations caused by laser instabil¬ 
ity, optical beam pointing drifts, and atomic motion. A 
small differential Stark shift also prevents time varying 
detuning, and consequently off-resonant state rotations, 
during the finite rise and fall times of the optical pulses. 
Because the experimental detuning is only several times 
larger then the 7 p 1 / 2 hyperfine splitting, the hyperfine 
structure must be included when calculating the differen¬ 
tial shift (see Appendix ©. State dep endent Stark shifts 
due to the 780 nm trapping light [la] are not accounted 
for since we turn off the traps for the few fi s duration of 
the Rydberg gates. 


B. Setting and Estimating parameters 

In order to calculate all of the shifts induced by the 
Rydberg excitation beams we need to have good mea¬ 
surements of the beam intensities at the atoms. A series 
of Ramsey and Rabi flopping experiments is used for this 
purpose as shown in Fig. [3] 


To find Si, two microwave 7 r/ 2 -pulses are applied with 
a variable length 459 nm pulse in between. The resulting 
Ramsey frequency is equal to the differential Stark shift 
of the qubit states induced by the 459 nm light. This is 
equal to the difference between the resonant Stark shifts 
on | 1 ) and | 0 ) which are separated by the hyperfine fre¬ 
quency w q , 

Adiff,g = A]q Aq X . (19) 

Measuring Adiff, fl , and using the dependence on the 
intensity of the 459 nm light we can infer S\ since Ai 
and w q are known. In principal the same method can 
be used to extract £2 but the differential shift on the 
ground state from the 1038 nm beam is small, about 
200 Hz. Instead we measure Dr the two-photon Rabi 
frequency for Rydberg excitation. Since Dr depends on 
£i,£ 2 , Ai, and the hyperfine structure of the intermedi¬ 
ate level, which is known, we can infer £ 2 . We emphasize 
that for our experimental parameters it is important to 
account for the hyperfine structure using the expressions 
given in Appendix [A] Although the full expressions only 
differ by about 10 % from the approximate expressions 
of Table |H the gate performance is very sensitive to the 
beam intensities at the 10% level, as can be seen in Fig. 
[7] below. 

As a consistency check we then measure the ground- 
Rydberg differential shift Adiff, g R with a ground-Rydberg 
Ramsey measurement. This shift depends on both £\ and 
£ 2 - Figure [3] shows measured data using the Rydberg 
825^2 state. The measured and inferred quantities are 
listed in Table El Using the matrix elements from Table 
IIIII we determine the field strengths £j which are used to 
calculate resonant and non resonant Stark shifts which 
in turn are used to infer Ai n f err ed.gR- Figure |3jc) shows 
the ground-Rydberg differential shift after choosing beam 
powers such that the shift is relatively small, less than 
100 kHz. At these low frequencies we are not able to 
measure the shift accurately, since we cannot hold the 
Rydberg atoms for extended periods. We therefore only 
give an estimated upper limit on the shift in the Table. 

Note the parameters above are not fine tuned to re¬ 
cover an ideal Cz as discussed in section lHDl We instead 
use a relative phase of <f >01 — 4 >00 between the ground state 
7 t /2 rotations to recover an entangling Cx,tc gate, as ex¬ 
plained in Sec. Ill El With the measured beam parame¬ 
ters, we use the equations of Sec. Ill Cl to calculate the 
expected two qubit operators 

( e -*o.i3 0 0 0 ^ 

0 e l0 - 88 0 0 

Cz,$= n7fi . ( 20 ) 

0 0 e -*0-76 q 

^ 0 0 0 —e* 0 24 ) 
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Ramsey gap time (ps) 




Rydberg pulse time (ps) Ramsey gap time (ps) 


FIG. 3. (color online) Experimental data used to determine optical beam powers. The figure shows measurements on the 
control site only and fitted frequencies are indicated in the plots, (a) ground state Ramsey oscillations with 459 nm laser 
used to extract Hi. (b) ground-Rydberg Rabi oscillations used to extract O 2 once Hi and Ai are known, (c) ground-Rydberg 
Ramsey experiment used to measure Adiff, g R. The straight line is an aid to the eye and is not a fit, but suggests a frequency 
< 0.1 MHz. The inset shows a faster ground-Rydberg Ramsey oscillation using other parameters. Beam powers extracted from 
these measurements are given in Table El 


measured values 

range over 2 days 

value used for fits 

value used for fits 

Cx eye diagram 


(MHz) 

control site (MHz) 

target site (MHz) 


Adiff ? g/27T 

0.80-0.95 

0.86 

0.81 


Or/2tv 

.63-.75 

0.67 

0.65 


Adiff ,gR /27T 

<(>01 ^ 000 

- 0.1 



1.05 (rad) 

inferred values 





^inferred, gR,/27T 

001 — 000 


0.088 

0.160 

1.01 (rad) 


TABLE II. Experimentally determined optical parameters. The first column give the range of measured values from multiple 
measurements over a two day time span. The next two columns give the values assumed for determining the optical beam 
powers. Using beam waists of W 459 = 3.0 /im, W 1038 = 3.7 fim the fitted powers were ^ 459,0 = 22. pW, Bisg.t = 21. /iW, 
-Pio 38 ,c = 1.9 mW, and Pio 38 ,t = 2.0 mW. The inferred ground-Rydberg differential Stark shift and Cx phase using these beam 
powers are given in the last two rows. 


and 


C x ^oi ~ </>oo) 


/ 0 e - * 1 - 7 0 0 \ 

e-iom 0 0 0 

0 0 e _l0 ' 76 0 

^ 0 0 0 e" 12 - 89 / 

( 21 ) 


All nonlisted elements in C x have magnitude < 8 x 
10~ 4 . We see that the gate phases differ appreciably 
from the standard values. Nonetheless we can still create 
entangled states as explained in Sec. Ill El The predicted 
Bell type state produced by this gate is 

w = c ^ nn m+e<- 2 *°°+^ + M\ll) 

^oo |00)+e l0 38 |H) 

V2 


which is a maximally entangled state. The largest ef¬ 
fect not accounted for in the predicted state is the fi¬ 


nite blockade strength of the two-qubit Rydberg inter¬ 
action. We use a Rydberg excitation Rabi frequency of 
Hr, = 2 -7T x .67 MHz. The Rydberg interaction for 82 si/ 2 
states at 7.6 /im separation is B ~ 27r x 23 MHz. We 
calculate a blockade leakage phase of 0 bl = ^ ~ 2.6 
deg., which is negligible compared to other error sources. 


C. Cx gate-experimental results 

It is clear that setting the correct phase of the second 
ground state pulse is crucial to the operation of the gate 
and creation of entangled states. This phase is found 
experimentally by varying the phase of the final ground 
state pulse. An example of this can be seen in Fig. [4] 
for control and target qubits that are two sites away so 
that their separation is 7.6 /im. Single atom data of the 
target atom which is cut on whether a control atom is 
present or not is shown. The blue curve shows the data 
when a control atom is present and therefore the Rydberg 
blockade occurs. While this is only single atom data, this 




























FIG. 4. (color online) Experimental Cx eye diagram on next 
nearest neighbor sites separated by 7.6 /im. The blue curve 
is from loading an atom in both control and target sites and 
represents a blockaded data set. The red curve is from data 
with no atom in the control site. The two curves have a n 
phase shift with respect to each other as expected. The phase 
at the minimum of the no blockade curve is q ioi — (f >oo which 
gives the Cx, tc gate. 



FIG. 6. (color online) Bell state population measurement and 
parity oscillation measurement with no loss correction using 
next nearest neighbor sites. The populations measured for 
the Bell state are |00) : 0.54 ± 0.06, |01) : 0.03 ± 0.02, 110) : 
0.05 ± 0.03, 111) : 0.38 ± 0.06. The fit to the parity curve 
gives |Ci| = 0.27 ± 0.02 and IC 2 I = 0.006. This results in an 
entanglement fidelity F = 0.73 ± 0.05 and F = 0.79 ± 0.05 
when corrected for atom loss. 

110), 111) the control qubit is Rydberg excited and must 
wait there while a Rydberg pulse is applied to the tar¬ 
get qubit, before returning to the ground state. Excess 
loss of population in this case is the largest contributor 
to gate error. We discuss possible reasons for this loss in 
Sec. ||V] below. 


state preparation 


|00> 

1.00 

.00 

.00 

.00 

|01> 

.00 

.97 

.00 

.03 

|10> 

.00 

.00 

.98 

.02 

I11> 

.00 

.01 

.01 

.97 


|00> 

|01> 

|10> 

|11> 


output state 


Cx populations 


|00> 

.07 

.88 

.00 

.06 

V 101 > 

.87 

.12 

.02 

.02 

$ 





■M 





Q-|io> 

.04 

.03 

.71 

.02 

|ii> 

.00 

.00 

.02 

.82 


| 00 > | 01 > | 10 > | 11 > 
output state 


FIG. 5. (color online) State preparation (left) and Cx gate 
population matrix (right) on next nearest neighbor sites, 
7.6 /rm apart, using parameters measured in Sec. IIII Bl The 
average statistical uncertainty of the data is ±0.008 for the 
state preparation and ±0.02 for the Cx gate. 


curve can be thought of as representing the 111) state and 
if we wish to run the C.y gate, we choose the phase where 
this curve is minimum so that 111) goes to 110). The red 
curve shows the data when a control atom is not present 
and no blockade occurs and can therefore be thought of 
as |01). 

To measure the Cx gate population matrix, each of 
the computational states are prepared, the Cx pulse se¬ 
quence is applied, and the results are measured as shown 
in Fig. [5] The overlap of the populations with an ideal 
gate is 0.82, without any corrections for atom loss. This 
improves on our previous result [fj of 0.74. Note that the 
sum of the output populations in a given row is not equal 
to unity. For the first two rows we get 1.01, 1.03. We 
attribute this slight excess to fluctuations in the load¬ 
ing and atom retention rates. The second two rows have 
populations sums of 0.80, 0.84. For these input states 


D. Entanglement Results 

We proceed to implement the U operator of Eq. ED 
to create an entangled two-qubit state. We select the 
gate phase of <foi — 4> 00 to implement the Cx gate. Data 
are taken on next nearest neighbor sites, separated by 
7.6 ^m. Parity oscillations are then performed to quan¬ 
tify the enaglement fidelity. The results of Fig. QT| give 
an entanglement fidelity, of F = 0.73 ± 0.05 without 
any loss correction. Adding in a retention correction for 
ground state loss during the gate sequence equal to 0.996 
and 0.993 for the respective sites does not change this 
result. Correcting for atom loss during the gate by re¬ 
normalizing each of the points on the parity curve results 
in a loss corrected value for Ci = 0.32 ± 0.03 and a post- 
selected entanglement fidelity of F = 0.79 ± 0.05. These 
results are the highest two-qubit entanglement fidelity, 
without post selection, reported to date using the Ryd¬ 
berg interaction. The previous best was 0.60 [111]. 


IV. DISCUSSION 

We have demonstrated improved entanglement fidelity 
using the Rydberg blockade interaction between two 
atomic qubits. Nevertheless the results obtained are still 
far from the 10~ 3 , or lower, errors that are expected to be 
needed for scalable quantum computing Detailed cal¬ 
culations do predict the feasibility of much higher fidelity. 
It is therefore important to understand the cause of the 
observed infidelity and to implement improved protocols. 

One issue is the technical challenge of stabilizing all 
experimental parameters. Our current experimental pro- 
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cedures for qubit state measurements involve pushing out 
/ = 4 atoms, and then detecting the presence of an atom, 
in order to infer the qubit state. This method allows us to 
measure the qubit state with high fidelity Q, but implies 
that a new atom has to be loaded half the time on aver¬ 
age. The need for atom reloading results in a relatively 
low data rate of 2 s -1 . While state dependent measure¬ 
ments have been performed without atom loss [ill, H5, 
to date this has only been demonstrated on one or two 
trapped atoms. In the multi-site array used here there is 
increased background noise from the trapping light and 
multiple atoms which makes lossless measurements more 
difficult. Achieving lossless detection in the multi-site ar¬ 
ray is a challenge we are working to solve by increasing 
the optical detection efficiency. 

The low data rate implies that entanglement experi¬ 
ments, including the requisite tuning of experimental pa¬ 
rameters, as in Fig. [3j require many hours to complete. 
As is shown in Fig. [7] the fidelity drops steeply as ex¬ 
perimental parameters are changed. In order to have a 
fidelity within 95% of the optimal value the beam powers 
should not differ by more than ±0.05 and the detuning 
should not change by more than ±0.2. It is not diffi¬ 
cult to maintain the detuning to the required precision, 
but holding beam intensities at the atoms to a few per¬ 
cent drift over many hours is challenging and needs to 
be improved on. Experimental intensity drifts are cur¬ 
rently up to the 10% level over the course of a day which 
contributes to the gate infidelity. 

In addition to these technical issues the measured data 
point to a dominant experimental error. Looking at the 
last two rows in the Cx population matrix in Fig. [5] we 
see that in cases where the control atom is Rydberg ex¬ 
cited there is about 20% loss of population. This error, 
together with imperfect state preparation, and drifts of 
parameters at the few percent level easily explains the im¬ 
perfect fidelity we obtain. Excess loss of Rydberg excited 
atoms has been the major contributor to gate infidelity 
not only in this study, but also in previous entanglement 
experiments We note that for input state 101) 

where the control atom is not excited, and the target 
atom experiences a 27 t Rydberg pulse, the loss is very 
low, at most a few percent as shown in Fig. [8] The dif¬ 
ference between excitation of control or target atoms is 
that the control atom has some gap time in between Ry¬ 
dberg 7r rotations where the excitation beam is switched 
off, moved to the target site, then moved back to the 
control site, and switched back on for the final Rydberg 
7 r rotation. This gap time is equal to ~ 3 /is. This is 
much shorter than the lifetime of the Rydberg state and 
only slightly longer than the 27 t Rydberg pulse applied to 
the target atom. We have verified in control experiments 
that a sequence of 7r Rydberg pulse, very short gap time 
of < 50 ns, 7r Rydberg pulse, also leads to excess Rydberg 
loss. 

The reasons for the additional loss when an atom is left 
in the Rydberg state for a short time are under investi¬ 
gation. Several possible explanations are worth consid¬ 


ering. The large polarizability of Rydberg atoms could 
result in mechanical forces from background electric field 
gradients pushing the atoms away. The experiments are 
performed in a pyrex vacuum cell with the atoms 1 cm 
away from the nearest walls. By way of spectroscopy on 
Rydberg states we have determined the background dc 
field to be 10 — 30 mV/cm. Although this measurement 
does not directly tell us about field gradients, assuming 
that the field varies by this amount over a length scale 
as short as 1 mm gives insufficient mechanical forces to 
explain the observed atom loss. Photoionization rates [dj 
are also too small to explain the observed loss. 

Another possibility is that what we observe is not 
loss of an atom, but coupling between the laser excited 
|82 Si/ 2 ,mj = —1/2} state and some other Rydberg state 
during the gap time. Since other Rydberg states are 
not brought back to the ground state in the second tt 
pulse, such coupling could result in leaving an atom in 
the Rydberg state, and the atom then being lost when 
the trap light is turned back on after the gate sequence. 
Several mechanisms could result in coupling between dif¬ 
ferent Rydberg levels. Although the nsi/ 2 states have 
no significant tensor polarizability, matrix elements to 
other states are very large, scaling as nreao. Higher or¬ 
der terms in the hypernolarizabilitv [2lj could then lead 
to state mixing. Improved control of background electric 
fields should serve to limit this possibility. On the other 
hand if a static background field was the reason for state 
mixing we might expect to see a similar loss in a single 
27 t Rydberg pulse, which we do not, as can be seen from 
Fig. E 

These considerations point to the likelihood that the 
observed loss is related to the turning on and off of the 
Rydberg pulses. The optical pulses are only approxima¬ 
tions to square pulses, as they have rise and fall times of 
about 20 ns. When there is a finite ground-Rydberg dif¬ 
ferential Stark shift there will be a time varying detuning 
at the rising and falling edges which may promote excita¬ 
tion of more than one Rydberg state. Figure^) shows a 
loss signal that does not increase monotonically with the 
gap time. Interestingly the minimal loss need not occur 
at the minimal gap time, and the population shows oscil¬ 
latory behavior. Although we have observed oscillations 
as a function of gap time, the details of the oscillations 
have not been repeatable. The oscillations may be due 
to time dependent interference of atomic states which 
are not eigenstates of the Rydberg excitation Hamilto¬ 
nian. Numerical integration of the Schrodinger equation 
for such pulse sequences does reveal imperfect Rydberg 
excitation and de-excitation, and this issue is also the 
subject of ongoing work. 

Irrespective of the correct explanation for the excess 
Rydberg loss it is apparent that tuning the ground- 
Rydberg differential ac Stark shift to a value that is small 
compared to the excitation Rabi frequency has reduced 
the amount of loss, and improved the gate fidelity com¬ 
pared to previous work. We anticipate further improve¬ 
ment in future experiments with better control of back- 
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FIG. 7. (color online) The entanglement fidelity as a function of the fractional change of the optical powers P459, P1038 and the 
detuning Ai from the optimal values. The center point of each plot uses the parameters found experimentally in Sec. IIII Bl 
and the color scale is the fidelity normalized by the value at the cener of each panel. 
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hyperfine structure of the intermediate level used for two- 
photon excitation. This is a good approximation when 
the detuning is very large compared to the width of the 
hyperfine structure. With our experimental parameters 
this approximation is only about 90% accurate. We give 
here relevant formulae that account for the hyperfine 
structure. 


FIG. 8. (color online) Ground state population after a 2tt Ry¬ 
dberg pulse a) and after a 7r— gap — n Rydberg pulse sequence 
with variable gap time b). 

ground electric and magnetic fields as well as the use of 
optimized pulse shapes. 
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We wish to couple atomic states |<?) —>■ |r) using two- 
photon excitation via intermediate level \p). When the 
hyperfine structure of the intermediate level is negligible 
the two-photon Rabi frequency is Or = ff iff 2/2 A with 
O12 the one photon Rabi frequencies and A = uq — 
LO P g the detuning of the field driving | g) —> \p). Here 
LUpg = ojp — ujg = (Up-U g )/h. This expression is accurate 
when |A| Ay where Ahf is the width of the hyperfine 
structure of the p level. 


Appendix A: Rydberg excitation with intermediate 
state hyperfine structure 

The expressions given in the main text for the ground- 
Rydberg Rabi frequency and AC Stark shifts neglect the 


We proceed to calculate the Rabi frequency for two 
photon excitation | g) -4 | p) -4 |r) where | g), |r) are spe¬ 
cific hyperfine states and | p) is shorthand for a manifold 
of hyperfine states. The two-photon matrix element be¬ 
tween ground and Rydberg hyperfine states is 


4/g — y fr = £ 1^2 e 2 (r, f r , m g + qi + q 2 \r q2 ^ | p, f p , m g + qi)(p, f p , m g + q 1 \r qi \g, f g , m g ) 

fp 

Jr fr r tfr,m g +q 1 +q 2 ^/p,m 9 +gi 
C Ijgfg C Ijpfp^fp,mg+q 1 ,l,q 2 ^f g ,mg,l,q 1 


(Al) 
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where 

V = £i£2e 2 {n r lrsjr\\r\\n p l p sjp)(n p l p sj p \\r\\n g lgsjg}, 
C " is a Clebsch-Gordan coefficient Q and 

4/ = {-i) 1+I+f+i ' Vv+i s j /,{j„ 


where Sf e c } = 



is a compact notation for the 6j 


symbol. Here e is the electronic charge, I is the nuclear 
spin, s = 1/2 is the electronic spin of an alkali atom, 
n g,p,r is the principal quantum number of the atomic 
state, /g lPl r is the total angular momentum, l g , p , r is the 
orbital angular momentum, and j g , p , r is the electronic 


angular momentum. The photon fields 1 and 2 have am¬ 
plitudes £d, (?2, polarization state q\, q 2 in a spherical ba¬ 
sis and m g is the projection of the ground hyperfine state 
angular momentum. 

We define a one photon detuning from the center of 
mass of the p state Ai = uq — u> pg and hyperfine shifts 
within the | p) manifold A f = uif — uj p = 27 r x ^[f p (f p + 
1) — /(/ + 1 ) — j p (j p + 1)] where A is the magnetic dipole 
hyperfine constant. We only consider states with j = 
1/2 and can therefore neglect higher order terms of the 
hyperfine interaction. 

Adding the contributions from the p levels the two- 
photon Rabi frequency between hyperfine states can be 
written as 

^fr,m g +qi+q 2 _ ,m g +qi+q 2 

fgi rn g fg: rn g 

where 


n 


n = 

fr,rn g +q 1 +q 2 
f g ,m g 


£i£2e 2 {n r l r sjr\\r\\n p l p sj p ){n p l p sj p \\r\\nglgsjg) 


I+jp 

= E 

f p —\i j p | 


2/i 2 Ai 


Jpfp Jrfr r ifr,mg+q 1 +q 2 /~if P ,m g -\-qi 


A! 


"Ijgfg Ijpfp fp,m g +qiP,q 2 f g , m g , l,9i ^Af ’ 

^ Jp 


r 


The Rydberg hyperfine states can be expanded in an 
uncoupled basis as 

\jrl;fr,m r )= Y C 'il™ r j im I \jrI\'rn j mi) 


rrij ,mj 


or 


\j r r, rrijmi) = Y <4™ 


j r mj Imi 


\jrI; fr^r)- 


fr,m r 


We are interested in the situation where we start in a 
single ground hyperfine state and therefore we replace 
m r by the laser excited value m g + qi + 92 ■ Then since 
= 2(/ r m T .|'H|/ g m g ), with H the electric dipole 
Hamiltonian for the two-photon transition, we can write 
the Rabi frequency coupling ground hyperfine and Ryd¬ 
berg fine structure states as 


C', = 


f r m r 




fr,m r 


j r rrij Imi fg,m 


fr,m 

= SJ E C ; 


fr,m g +qi+q 2 f.f r ,m g +qi+q 2 
jrmjlmr ^ l fg,m g 


= QQ Jr ’ mj 

~ UU fg,m g 


(A2) 


with to/ = mf g +ql + q2 — rrij r . Here we have introduced 
an effective angular factor 


O 


n 


E ^fr,m g + l 

°jr ,rn,j ,1 ,m g +q 1 +q 2 -m j i 11 f g ,m g 
fr 


fr,rn g +q 1 +q 2 


(A3) 


We use tildes to denote angular factors coupling hyperfine 
states to hyperfine states and overbars to denote factors 
coupling hyperfine states to fine structure Zeeman states. 
It should be emphasized that a description of the ground- 
Rydberg coupling in terms of a Rydberg fine structure 
state is only valid when the hyperfine interaction in the 
Rydberg state is negligible. Cesium ns states have a rel¬ 
atively large hyperfine splitting so that in order to ensure 
coupling to a single Rydberg state we apply a bias mag¬ 
netic field along z to decouple the hyperfine interaction. 
This implies a modification to Eq. (IA3I) which we will 
make explicit in Appendix [Bl 

We will also need the one-photon Rabi frequencies 


^fp 

f g 


77 Jfv 
^9P^f g ,m fg ,q 1 

77 Jpfp s~lfp’ rn fg-^Q 1 


^IP 

^3r,m jr -q 2 


^ r P"j r ,m jr ,-q 2 


— £}. r P Y] 

fr 


Jpfp /^ifp, m f r 92 r<fr-:<rif r 
C Ijrfr fr,m fr ,l,-q 2 ^ / j r m jr Imi 


where m/ r = to/ s + qi + <72 , to/ = m/ r — m,j r and 

,—, _ £\&(ji p l p sj p \ |r| \n g l g sjg) 

“ 9P = h ’ 

,—, £ 2 &(n r l r sj r | |/j \n p l p sj p ') 

= li ' 

With these definitions the ground and Rydberg state 
Stark shifts under conditions of two-photon resonance 









12 


are 


(a f p f 

' aC ’ g ^ 9P 4(A X — A/ p ) 

(h* Y 

_ V Or,m jrt -q 2 J 

lac ’ r_a ^4- 4 (Ai-A /p ) 


(A4) 

(A5) 


The differential AC Stark shift is A ac = A ac>r — A aCjg . 

We have written the Rydberg fine structure state AC 
Stark shift as a sum over the contributions from hyper- 
fine states of the p level. When the Rydberg hyperfine 
coupling is not negligible we should instead use the hy¬ 
perfine resolved AC Stark shifts which are 


„ 2 y- (“A,™/,.,-92) 
“ rp 4' 4 (A x -A fv ) 

JP 


Finally, it is also important to know the probability 
of photon scattering from the p level during a ground to 
Rydberg tt pulse. The time for the pulse is t v = 
and the number of scattered photons is (see Sec. IV.B in 

0 ) 


r y 

*T _ 7 pin —2 V /s- m fl-9i/ 

2 a «'4-2(A 1 -A /p )2 



afp ) 

3r,m j -q 2 J 


21 


2(A 1 -A /p ) 2 


Here = 1 /t p is the radiative decay rate from the p 
level and the prefactor of 1/2 acounts for half the in¬ 
tegrated population being in the ground and Rydberg 
levels during the 7r pulse. For coherent qubit control we 
choose parameters such that N -C 1 and interpret N as 
the probability P se to scatter a photon. 

The above general expressions for Rabi frequency, AC 
Stark shift, and spontaneous emission probability can be 
applied to any desired set of atomic levels. Expressions 
valid in the limit of detuning large compared to the hy¬ 
perfine structure width can be found by simply setting 
A f -A 0. 

JP 


Appendix B: Excitation of Cs nsi /2 states via 7pi/2 


For the experiments reported here we use Rydberg 
level 82s 1/2 which has a hyperfine splitting between the 
f r = 3,4 components of about Whf.r = 27r x 110. kHz. At 
zero magnetic field the |6s!/ 2 ,4,0) ground state couples 
to both hyperfine levels with equal and opposite Rabi fre¬ 
quencies. The non-zero hyperfine splitting implies that 
for any laser tuning there is some off-resonant excitation 
which leads to gate errors. 


parameter 

value 

Ref. 

nr 

«g,459 

-11.6 x lO " 24 (cm 3 ) 

a) 

nr 

a g,1038 

189. x 10 “ 24 (cm 3 ) 

a) 

(7pi/ 2 |k||6si /2 ) 

-0.276 a 0 

M 

A 7 P 1/2 

94.35 (MHz) 

m 

A 3,7p 1/2 /(27r) 

-212.3 (MHz) 


A 4,7p 1/2 /(27r) 

165.1 (MHz) 


r lPl /2 

0.155 (ps) 

[26] 

a 82s 1/2 ,459 

-15. x 10 “ 24 (cm 3 ) 

Eq. © 

nr 

Q! 82s 1/ 2,1038 

-77. x 10 “ 24 (cm 3 ) 

Eq. © 

(ns 1/2 \\r\\7pi/2) 

0 

b) 

AnSi/ 2 

FtHf ( MHz ) 

c) 

^82 S 1/2 

203. (ps) 

m 


TABLE III. Physical parameters for Rydberg excitation of 
Cs via the 7 pi/2 level. From top to bottom the table sections 
give ground state parameters, 7 pi/2 parameters, and Rydberg 
level parameters. Reduced matrix elements are given in terms 
of the Bohr radius ao- 

a) The ground state nonresonant polarizabilities are calcu¬ 
lated using a sum over states method, excluding the 7 pi /2 
level in the case of 0 ^( 459 . 

b) The n dependence is a fit to values calculated using quan¬ 
tum defect wave functions as described in [ 2 §], 

c) The n dependence is an approximation based on values 
reported in [l5|, |2!|■ 


Here we give explicit expressions for excitation of Cs 
atoms using 6s!/ 2 —> 7pi/ 2 -A ns 4 / 2 . The two fields 
are er_|_,<r_ polarized so that the ground \Gsi/ 2 ,f g = 
4:,rrif = 0) state is coupled to the Rydberg hyperfine 
states |nsi/ 2 ,3, 0), |ns!/ 2 ,4,0). The hyperfine-hyperfine 
angular factors defined in Appendix [A] are 


~3 

■"4,0,1 


53 

“3,0,1 


~3 

“4,0,D 


,4 _ VW3 A4 

“4,0,1 — 4 ) “3,0,1 


^4 

'“4,0,1) 


Q 4 >° — 
“4,0 ~ 


1 


Ar 


16 Ai - A 3 


5 q3,0 _ 

48 Ai — A 4 ’ 4 ’° 


-fi 


4,0 

4,0' 


To correct for this we apply a small bias magnetic field 
B z along the z axis of B z = 0.15 mT. In the presence of 
the magnetic field the energies of the Rydberg |3,0), |4,0) 
states move apart and the coupled eigenstates can be 
written as 


182s 4 / 2 , rrij = 1/2) = 


|82 s 1 / 2 ,m j = -1/2) = 


- (l - VlTx 2 ) |3,0) + x|4,0) 

_ n -| 1/2 : 

x 2 + (1 — VTTx 2 ) 

(i + VTT^0 |3,0) — x|4,0) 


c 2 + (l + a/1 + x 2 ) 


11/2 


with x = pBgjBz/hur,hf where ps is the Bohr magneton 
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and gj ~ 2. Here we have neglected the small correc¬ 
tion due to the nuclear g factor. At B z = 0.15 mT we 
find x = 38.2 and the matrix element from the ground 
state to |82s 1 / 2 ,m J - = —1/2) is within 0.0001 of the 
asymptotic value. At the same time the coupling to 
|82 si/ 2 ,TOj = 1/2) is suppressed to about 0.02 times 
the coupling to |4,0). The small matrix element to¬ 
gether with the approximately 4.2 MHz splitting between 
the rrij = ±1/2 states lets us tune to resonance with 
rrij = —1/2 and safely neglect the coupling to rrij = 1/2. 
Measurements indicate a residual coupling to rrij = 1/2 
that is slightly larger than 0.02, which may be attributed 
to polarization errors of the Rydberg excitation beams. 

To summarize, with the magnetic field applied we cou¬ 
ple to a Rydberg fine structure state with matrix ele¬ 
ments given by Eqs. (IA2IA3IA5I) . Using the expressions 
in the previous section we find the “reduced” one- and 
two-photon Rabi frequencies 

= 87570. x J-P^wP 2 ,mw) 7 - ( MHz ), 

w pl/2 

= 2446. x 1,mW (MHz), 

27T W l, M m 

r-n pi / 1 2 

^ = 71600. x 2 ' mW (MHz). 

2?r n r W 2 l(lm 


where vi^ghz = Ai/27t (in GHz), Pj , m w is a beam power 
in mW, and w\( 2 ),pm are beam waists (1/e 2 intensity 
radii) in pm. The two-photon angular factor coupling 
ground hyperfine to Rydberg fine structure state can be 
written compactly as 


5A 3i7Pi/2 3A 4 ,7p 1/2 

_ 8 A 1 

)0-^) 


and ST 3 /g 2, 1 ^ 2 = — H^q 2, . Finally, the full expressions 

for the ground-Rydberg Rabi frequency, AC Stark shifts, 
and spontaneous emission probability from the 7p\/2 level 
are 


O 7 • 7 

“ 4,0 


3x/2 


; ’ 7 P 1/2 


Ai 


n 


1 / 2 ,- 1/2 
4,0 


27T 


= -20600.- 


( p i, 


mWb,mW 


)l/2 


1 - 


5 A 3; 
8Ai 


3A 4j 

8 A 1 


Aac 

27T 


nr /2 Wl, M mW2, M m^l,GHz (l - (l - 

^ac,r ^ac,g 


27r 


(i 


= I 160200. ^ 2,1 ” w - 93.47 Pl ’ mW 


1 




^2,/ira 


»l, M m 


= 


0.43 ( P2 ; mW ) 1 2 + 2.5 X 10 - 4 

W 2 , M m 


l — 2 ( 5A3p 
1 ^ 8Ai 


8Ai J 


+ 


5A, 


"8A? 


3A 4 

*P 

8A'{ 


(l-^)(l-^)(l 


5A 3 , 

8Ai 


1 


t'l.GHz \ 1 — 


Aa i 

A, 


SA^V 
8Ai J 


(MHz) 


5/3 


1 ^ 4 P 

1 Ai 


v P2,mW 


(MHz), 


1 


»2,(i m 
Wl.jim ^l,GHz 


The ground-Rydberg Rabi frequency, resonant ground 
Stark shift, and total ground-Rydberg Stark shifts are 
shown versus detuning using experimental parameters in 
Fig. [9] Comparing the center and right panels demon¬ 


strates that the non-resonant Stark shifts must be ac¬ 
counted for when determining the parameters to be used 
to cancel the ground-Rydberg differential shift. 
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FIG. 9. (color online) Calculated ground-Rabi frequency Hr (left), resonant Stark shift of | 6 si/ 2 , f = 4) due to the 459 nm light 
(center), and total ground-Rydberg differential Stark shift A g R including resonant and non-resonant contributions (right), 
for two-photon excitation of Cs |82 s,rrij = — 1 / 2 ) starting in |4, 0 ) using parameters from Table ITT1 for the control site. Solid 
blue curves account for the 7 pi /2 hyperfine structure, while the dashed yellow curves do not include hyperfine corrections. The 
spontaneous emission plot only accounts for scattering from 7 pi /2 as the Rydberg state decay is negligible for the parameters 
used. The gate experiments used v\ — 0.83 GHz. 
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